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COMPARATIVE STUDY OF FLARE CONTROL LAWS 


By 

Arun A. Nadkarni* 

INTRODUCTION 

This report presents the development o£ a digital, 3-D, automatic 
control law designed to achieve an optimal transition of a B-737 aircraft 
between glide slope conditions and the desired final touchdown condition. 
The digital control law developed here is a time-invariant, state-estimate 
feedback law, and the design is capable of using the Microwave Landing 
System (MLS) under development by the Federal Aviation Agency (FAA). The 
study of a curved flight path leading to a steep final approach and touch- 
down under low visibility conditions is part of the Terminal Configured 
Vehicle (TCV) program, sponsored by NASA/Langley Research Center (LaRC). 
The goals of this program include the reduction of aircraft noise in 
communities surroundings airports, Che reduction of fuel consumption, the 
reduction of the effects of adverse weather conditions on aircraft opera- 
tions, and the efficient use of airspace in congested terminal areas. A 
specific objective which supports these goals is the development of the 
capability to perform automatic flares from steep glide slopes to precise 
touchdown locations. 

The major reason for the use of steep glide slopes is the resultant 
noise reduction in comparison with the currently used 2.3* to 3* glide 
slopes with the Instrument Landing System (ILS). The steeper glide slope 
reduces the noise levels perceived on an identical segment of the ground 
for two reasons: first, at equal distances from the touchdown point, the 

aircraft flying a steeper, say a six-degree glide slope, is at about twice 
the altitude compared to that flying a three-degree glide slope. This 
difference in altitude causes a considerable reduction in the noise 
level perceived on the ground even if the two sources generate identical 

'^'Visiting Assistant Processor, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508, 
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noise: levels* Secondi en aircraft flying a steeper glide slope requires a 
lower! thrust setting, and this causes further reduction In the noise level 
perceived on the ground* The reduction in thrust setting has the added 
advan^ges of reducing the fuel consuoptlon during the final phase of the 
flight path* The ability to fly varying glide slopes may also provide an 
effective method to avoid encountering vortices generated by larger 
aircraft* This versatility toay result in a more efficient use of the 
airspace* i 

This Is not without Its attendant disadvantages* The use of steeper 
glide' slopes for noise reduction requires that the aircraft be flown In a 

i 

hlgh-drag, low-power-settlng configuration* In addition to this situation, 
the higher sink rates associated with these paths allow pilots considerably 
less reaction time to recognize an emergency and to take appropriate 

corrective actions In the event of atmospheric disturbances (e.g. , gusts, 

> 

wind shears, etc*), system or sensor failures, etc* 

t i 

phere is also a need for reducing touchdovm dispersions in the 
presence of varying flight conditions encountered during the flare portion 
of the landing. Ihe reduction In the touchdown dispersions greatly 
facilitates high-speed rollout, which significantly Increases the traffic- 

; j 

handling capacity of the terminal* It is obvious that developing 
capability to perform automatic flare maneuvers will accomplish many of 
these goals* i 

'In the next section, "Description of the System Model,” the system 

i i 

equations of motion of the aircraft are presented. These linear equations 
represent the perturbed motion of the aircraft from the nominal glide slope 
trajectory. A method of incorporating the spatial, low-level wind shears 

into I these perturbed equations of motion is indicated* It Is shown that 

i i 

the system equations Chen assume the familiar form of the linear regulator 
problem, acted upon by a constant disturbance. Under "Derivation of the 
Control Law," a design procedure to compute a digital, time-invariant, 

optimal control law for the discrete regulator problem acted upon by a 

i i 

constant disturbance is indicated* This is followed by a section 
desctlblng implementation of the cont.rol law. Under "Results," performance 
curves are presei^ted to show the capability of this digital, state-feedback 
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controller to perform the optimal flare maneuvers in the presence of 
various wind shear conditions Indicated. In the final section the 
conclusions derived from the study of the automated flare maneuvers are 
listed along with scope for further work. 

DESCRIPTION OF THE SYSTEM MODEL 
Introduction 

The development of the mathematical model In this study follows 
closely the development of a similar model described in detail in 
references 1 to 6. The complete derivation of the system equations Is 
described In these references; however, a brief outline of the derivation 
is given below for the sake of completeness. 

Aircraft Dynamics with Wind Disturbances 

This study Is concerned only with the final phase of flight, viz the 
flare. Thus, the aircraft is approaching the runway on a certain initial 
glldepath. The aircraft Is aligned with the runway, has a zero or at most 
a very small yaw angle with respect to the runway as well as a zero bank 
(or roll) angle, except in the case of a significant crosswind. Tlterefore, 
all the lateral dynamics are neglected during the analysis and only 
longitudinal dynamics are considered. 

With these assumptions and assuming small perturbations about the 
nominal path, the nonlinear equations of motion of the aircraft can be 
linearized using well-known methods. The complete equations of motion and 
the linearization procedure are outlined in references 5 to 7. These 
nonlinear equations are derived assuming (1) a flat earth, (2) an earth- 
fixed frame of reference, (3) a rigid aircraft, and (4) that second-order 
i.erms are neglected. These equations of motion are coupled. However, for 
a steady-state flight condition, the equations can be decoupled into two 
groups, the longitudinal equations and the lateral equations of motion. As 
already indicated, only longitudinal equations of motion are dealt with in 
this study. 


3 



The decoupled, nonlinear longitudinal equations of motion are then 
linearized about the nominal trajectory (!••, the steady-state flight of 
“3*, -6* glide slope, etc. ) to obtain the linear perturbation equations In 
the state variable form. The equations are expressed In a stability axes 
coordinate frame attached to the aircraft at the center of mass (fig. 1). 

The final linearized, longitudinal equations of motion, Including wind 
disturbances, assume the following standard form (refs. 1-A)*, 

X “ Ax + Bu + Dm (la) 


where 


X “ (C u' o q — n — 2. — g — 2, 6T 6th 6s 6e)^ 

(lb) 



u “ (5e 6s 

(Ic) 


and 

w •• (u^ Oy q„)^ (Id) 

where the states (x's) and the controls (u's) are 

6 ■ perturbation In pitch angle 

u' perturbation In velocity along x^ (stability) axis 

normalized 

o ■ perturbation In angle of attack 

q ■■ perturbation In the pitch rate 



perturbation In the horizontal position of the aircraft 
(normalized). In Inertial frame 



perturbation In the vertical position of the aircraft 
(normalized). In Inertial frame 


6T ■ perturbation In the thrust 


dth ■ perturbation in the throttle position 


5s ■ perturbation In stabiliztr position 


A 




fie “ perturbation in the elevator position 

fie *• perturbation in the elevator rate 

fiS ■> perturbation in the stabilizer rate 

fitli ■ perturbation in the throttle rate 

i , 

5sp ■ perturbation in spoiler position 

The subscript w indicates the perturbation in the variable due to 

wincj disturbances. Note that, out of the 10 variables in the state vector 

I 

X, the first 4 variables are sufficient to describe the longitudinal 

X - X 

perturbed motion' of the aircraft. The next two variables viz — 


z - r 

, ' o 


rrei defined for designing a 4-D control law to minimize the 


deviations from the nomina?. path, if desired. The perturbation in the 

I . i 

thrust IS defined as a state variable in order to model the thrust dynamics 
taking the "spool up" time of the engine into account, at least linearly. 

i : 

The ’last three variables, the perturbations in the position of the throt- 
tlei the stabilizer, and the elevator, were added as a result of the design 
decision to command the rates of these controls. 

I i 

{ The position bf the spoiler was included as a control variable in case 

it is decided to study the effects of direct lift control in future studies. 

In the present work, however, the spoiler was made inoperative during the 

. I j 

simvilation runs . ‘ 

Wind Model 

In order to: complete and simulate the system model given by equation 
(l),i the wind perturbation vector w muse be specified. The components 

I ^ I 

of this vector consist of u' , the normalized wind velocity in the +x 

direction; a , the perturbation in the angle of attack due to the wind; 

! ; 

and! q , the perturbation in the pitch rate of the aircraft due to wind. 

I i 

These wind variables may be logically modeled as the sum of a gust corapo- 

! i 

nent with zero mean value and a steady component. 

j The gust components are modeled using the well-known Dryden spectrum 
(re^s. 5-7). This method consists of using spectral factorization methods 

to obtain a dynamical system which generates a random process having the 

i j 

specified power spectral density when driven by white noise. Because of 


i 

i 


the linearity of the system, the three gust components can be treated 
individually: thus, only appropriate components are included in the 
longitudinal equations (eqs. (la) - (Id)). The detailed derivation of the 
gust components can be found in references 3 to 6 and will not be given 
here. 

The Dryden spectra describe the statistical behavior of the wind gust 
velocities in the aircraft body-fixed coordinates, and the gust components 
can be expressed in the following form (refs. 3-7): 


W - A W + Be, e, 
g ww g 1 


where 


g gb gb ^gb gb 


( 2 ) 


These four gust components constitute the four components (x^ to 
of the state vector. The elements of and can be 
obtained from reference 3. 

The steady-state components of the wind are simpler to model since 
they do not involve spectral factorization. These components can be 
modeled as the output of a first order deterministic plant, corrupted by a 
white noise. This can be done in different ways (refs. 3 and 4). 

To permit modeling of spatial shears, a new method was proposed (refs. 
8-10). The model also makes the whole system controllable, and the 
feedback gain matrix can be computed as is done in the usual manner for the 
regulator problems without the need for splitting the matrix Riccati 
equation in two. 

Let the aircraft be coming down on a nominal glide slope in a steady 

flight condition. Also, for simplicity, let the aircraft be flying in a 

disturbance-free atmosphere {t • s, s <o}. (It is noted that the method 

can be trivially extended to the situation when the aircraft is flying in a 

constant wind, for {t “ s, .<! <o}.) Now, at t “ 0 let the aircraft encounter 

a step wind, with a component U in the ^x direction and a component 

'^o ® 

W in the +Z direct ton. Also, let the subsequent wind field be 
o 




described by a linear shear profile with a shear rate of X kn/30. 46 m 

(100 ft) In the horizontal wind velocity, and Z kn/30.48 m (100 ft) in the 

vertical wind velocity. In the present notation, U >0 represents a 

o 

tailwind, W >0 represents a downdraft, X > 0 represents a linear 
o 

increase in tailwind (with decreasing height), and Z > 0 represents a 

linear increase in downdraft (with decreasing height) as seen by an earth- 

fixed observer. Then the rate of change of the perturbation wind velocity 

in the horizontal direction (after normalizing with ) is given by (ref« 

o 

10 ) ; 

* X 

u' ■ - TTrw (cos y ‘S + sin y «u’ - cos y «a) (3) 

W lUU 0 0 0 

Similarly, the rate of change of the normalized perturbation wind 
velocity in the vertical direction is given by 

• z 

w' «• - _irr (cos y *6 + sln y ‘u’ - cos y * 0 ) (4) 

w . lUU o o o 

Note that equations (3) and (4) yield the perturbation components of 
the wind at the wing. However, due to inertia effects, the flow field over 
the wing does not sense these changes until it travels a few chord lengths 
(ref. 5). These inertia effects can be modeled by a first-order lag term 
with an appropriate time constant. Thus, the steady- state components of 
the perturbation wind velocity in a shear field can be modeled as 

*15 " *1*15 **1*17 

*16 “ *1*16 **1*18 

Jt.T ■ “ (cos y • 6 + sin y • u* - cos y_ • a) 

17 100 00 o 

*18 “ ■ TWT ^o * ® y^ . u« - cos y^ • a) (5) 


where 








1 

I _ ; 

X '■ ** horizontal component of actual perturbation wind veljcity 
li felt by the wing, u' 


X i ■ vertical component of actual perturbation wind velocity 
ijS felt by the wing, w' 


X I ■ horizontal component of perturbation wind velocity at the 
17 wing, u' i , and 


xi “ vertical component of perturbation wind velocity at the 
18 wing, w'; 


The 4 steady-atate ^variables can now be augmented with the 10 aircraft 
statesj (eqa. (la) - (Id)) and A gust components to yield the complete state 
variable model of the aircraft motion in the presence of atmospheric 

'I . . 

disturbance. j 

i ; 

The system model developed above is completely controllable in 
addition to. being more realistic than the previous attempts to model the 

i 

atmospheric disturbance effects. The optimal, constant feedback gain 
matrix for this system can now be computed in the usual way as for a 
regulator problem,,' with the two nominal shear components acting as constant 
disturbances on the system. The method to compute the gain matrix is 
described in detail under "Derivation of the Control haw." 

i 

The Basic Augmented System: Aircraft in c. 

Wind Disturbance Field 

i(ith the derivation of the wind model for steady wind velocity 
components now complete, it is possible to augment the wind model (eq. 
(5))with the aircraft equations of motion in presence of wind disturba.ices 
(eqs.j (la) - (ld)j to yield the state equations of the basic augmented 


system. j 

iThe final form of the wind model assumes the following form: 


8 


( 6 ) 


W - A W + Bj., C, 
g ww g Cl 1 


8 WX 


A __ X «d 


(7) 


where C's are white noise processes to account for the ^mknown distur- 
bances and the wind vectors are 


W “ (a u ® u V 
g gb gb ^gb gb 


W 


(u* w' u* w' )'^ 
WWW w 
a a 


Here the subscripts g and s refer to gust and steady components, 
respectively. The elements of the matrices A^, , and 

can be obtained from references 3 and 4. The elements of the A^ 
and B^ can be obtained from equation (5) and references 9 and 10. 

Defining a composite wind vector 


W A (W^ (8) 

“ g s 

the wind vector w defined in equation (1) can be expressed as 



where (c „ ] is an appropriate transformation matrix (ref. 3). 

Wg w8 

The steady winds, the wi.nd shears, and the gusts can now be included 
in the system equations (la) to (Id). 

The complete system equations, with the inclusion of the wind model, 
can now be expressed in the stand^'^rd state variable form as 








>:“Ax + Bu + B W. + Bp5 
w d ^ 


, DERIVATION OF THE CONTROL LAW 


It ia now possible to undertake the design of the optirosl control law 

I 

for the above problem by invoking the separation theorem in the usual 
manner. The general approach to study the flare performance (refs. 1, 2, 

8) is briefly as follows. The longitudinal equations of motion of the 
aircraft (eq. (10)), which are perturbations about the nominal glide slope 
trajectory, are discre^zed (ref. 11). The constant gain, state feedback 
optimal control law, designed in a manner described below, is incorporated 
into the system equations, and the system is expressed as a deterministic 
closed-loop system. The difference between the initial glide slope and the 
desired touchdown conditions is supplied as initial condition x^, and 
the time response of the deterministic closed-loop system is simulated on a 
digital computer (refs. 8, 9, 12). 

The qomplete system equations for the augmented system can be 


expressed 


in discretized form as (refs. 8-10); 


The ^ptimal time-invariant control law for this system can now be 
computed jjsing the metjhod described in references 13 to 15. The control 
law can be expressed as 


1 - * H, , X, + H, Ju, 


f * W) 




Substituting for u. as 
•'k 

“k+l " “k .. .. . . . 

Uj « ■ ■, AT ■ txmc interval 

Vl - «ll\ " * '"12 " ' «“k • "l2” “k.O ■ “O 

The initial condition can be obtained from the expression for the 
minimum performance index (ref. 16) as 

—1 T 

uq ■ (P22 P12 *0 W) 

vhere (14) 

?n Pl2 

P - T 

P12 P22 


is the steady-state solution of the Riccati equation. This matrix is 
already obtained during computation of the constant, state feedback gain 
matrix above. 

Substituting the above control law, the original system equations can 
now be solved by simulation on a digital computer as usual, 


IMPLEMENTATION OF THE CONTROL UW 


Introduction 

In this section, two methods of implementing the digital, state feed- 
back, optimal control law derived in the previous section are discussed 
briefly. Their relative advantages and disadvantages are discussed under 
"Conclusions." 


Forced-Regulator Method 

The dichotomy of this method of implementation can be described as 
follows. Hy running a few’ trial simulations, the values of the rates of 
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throttle and elevator are found which will flare the airplane from the 
nominal glide slope (**3* in the present case) to an acceptable touchdown 
condition in the absence of any winds, either gusts or shears. Ttte values 
of these ramps on the throttle and elevator are now stored and forced on 
the airplane (in an open-loop manner) even wlien it is flaring through a 
given wind shear field. However, the perturbations in the state variables 
from the baseline no-winds trajectory due to the effect of the wind shears 
are now used to close the loop and generate an additional control using Che 
control law derived in the previous section. The total control therefore 
consists of two parts: (1) the open-loop ramps (i.e. values of the rates of 

throttle and elevator) derived for the no-wind flare performance and (2) 

the state feedback, optimal control law designed to drive the deviations 
from the baseline no-wind trajectory to zero. The two parts are added 
algebraically to yield the total control required. The particular 
structure of the controller leads to the name "Forced-Regulator Method." 

Regulator Method 

The concept of this method is very simple. The differences between 
the terminal (touchdown) conditions and the initial glide slope are 
supplied as the initial conditions for Che optimal, state feedback control 
law derived in Che previous section. The terminal conditions used could be 
either those on an equilibrium trajectory with sink rate of 0.61 m/j (2 

ft/s) or any desired touchdown conditions. The latter case was extensively 

studied in previous works of the author (refs. 8-10, 12). 

A few runs were simulated using the terminal conditions on an 
equilibrium trajectory with a sink rate of 0.61 ro/s (2 ft/s) to derive the 
initial conditions for generating Che closed-loop control law, but Che 
res ilts were found to be very unsatisfactory, and so are not presented 
here. The disadvantages are discussed in detail under "Conclusions." 

RESULTS 

In this section the performance of the discrete, optimal control law 
implemented using the Forced-Regulator Method described in the previous 
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section is evaluated* 'Jlie plant dynamics used are chose of the TCV Boeing- 

j 

737 research aircraft at NASA/LaRC. The aircraft Is assumed to be on a 3* 
gll!de slope when the flare maneuver Is initiated at a specified height. 
Simultaneously, the aircraft enters one of Che example shear profiles 
tabulated In cable 1. 

{ Note Chat, even chough an optimal feedback gain matrix was computed 
for! each case, for purposes of simplicity In onboard Implementation of the 
corjtrol law. It was decided to compute only one gain setting for the shear 

i 

profile A, and Co use Che same setting for all the cases simulated. This 

procedure would eliminate Che need for changing the gain setting onboard 

i 

Che: aircraft, each time a different shear profile Is encountered by Che 
aircraft. The time response curves presented Indicate the performance of 
the; Forced-Regulator Method of implementation of the control law described 
earlier. I 

I After a few Initial trial simulations. It was found that. In Che 
absence of any winds, ramiJing back the elevator at the rate of 0.7835 deg/s 
and the throttle at the rate of 4.193 deg/s yielded a very satisfactory 
flare performance (fig. 2). Touchdown occurred at 6.2 s at a horizontal 
distance of 393.8 m (1,292 ft). The sink rate at touchdown was 0.67 m/s 

i 

(2.186 ft/s) and the pitch attitude was 2.662°. The value of the thrust 

j : 

was reduced to 10,822 N (2,433 lb). It was noted, however, that the value 
of ;the sink rate at touchdown was very sensitive to the comand rate of the 
ele|vaCor. I 

With the flare trajectory thus obtained by forcing the ramps on the 

I 

elevator and the throttle as a nominal or baseline trajectory, a wind shear 
field of profile A (cable 1) was forced on the airplane. The state feed- 

( I 

ba(jk optimal control law was activated to generate an additional control to 
null the deviations of Che aircraft from the nominal trajectory (fig. 2). 
Figures 3a to 3c show the flare performance aa the aircraft flew through 
Chej shear profile A. The touchdown occurred at 6. 1 s at a horizontal 
distance of 38/j7 m (1,272 ft) from the Initiation of the flare maneuver, 

at la pitch attitude of 4. 387°. The sink rate at the touchdown was 0.66 m/s 

1 ! 

(2.;163 ft/s) an4 the flare was initiated at 15.54 m (51 ft) above ground 

i I 

level. 


I 

I 

I 
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Table !• Inertial Wind Profiles Simulated' 


Profile 

A 

B 

C 

D 

E 

F 

G 

H 


Hind Parameters 

Tailwind Increasing at 10 kn/30.48 m (100 ft) 
Downdraft increasing at 2 kn/30.48 m (100 ft) 

Headwind increasing at 10 kn/30.48 m (100 ft) 
Updraft increasing at 2 kn/30.48 m (100' ft) 

Headwind increasing at 10 kn/30.48 m (100 ft) 
Downdraft increasing at 2 kn/30.48 m (100 ft) 

Headwind increasing at 3 kn/30.48 m (100 ft) 
Updraft Increasing at 2 kn/30.48 m (100 ft) 

Tailwind increasing at 5 kn/30.48 m (100 ft) 
Dof^ndraft increasing at 2 kn/30.48 m (100 ft) 

Tailwind increasing at 15 kn/30.48 m (100 ft) 
Downdraft increasing at 2 kn/30.48 m (100 ft) 

Headwind increasing at 15 kn/30.48 m (100 ft) 
Updraft increasing at 2 kn/30.48 m (100 ft) 

Hecdind increasing at 15 kn/30.48 m (100 ft) 
Updraft increasing at 5 kn/30.48 m (100 ft) 
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Figure 3b shows that the engines were throttled down to reduce the 
thrust at touchdown to 10,924 N (2,456 lb). The elevator assumed a nega- 
tive value to provide the necessary pltch-up attitude. 

Figure 3c shows the velocity components of the wind at the wing and 
those felt by the wing. 

Figures 4a to 4c demonstrate the flare maneuver performed by the air- 
craft entering Into a wind profile B. It Is again noted that, when gener- 
ating the closed-loop control law for the shear profiles B through H, the 
optimal gain matrix computed for the shear profile A Is used In all the 
cases, even though. In each case, the optimal gain matrix was computed and 
found to be significantly different for each case. 

It Is clear that the reverse of the shear In both the horizontal and 
the vertical wind velocity components had a significant effect on the 
elevator history. The touchdown occurred at 6.0 s at a distance of 381.3 m 
(1^251 ft) from the Initiation of the flare. The pitch attitude at 
touchdown was 1.404®, the sink rate was 0.607 m/s (1.992 ft/s), and the 
thrust was 10,644 N (2,393 lb). It Is found that the flare would have to 
be Initiated at 14.66 m (48.1 ft). 

Figures 5a to 5c Indicate the flare performance when the aircraft 
entered the wind profile C. The touchdown occurred at 6.1 s at a distance 
of 386.8 m (1,269 ft). At touchdown the altoraft had a pitch attitude of 
2.513®, a sink rate of 0.658 m/s (2.162 ft/s) and a thrust value of 12,557 
N (2,823 lb). The flare was initiated from an altitude of 15.21 m (49.9 
ft). 

Figures 6a to 6c Indicate the flare performance when the aircraft 
encountered the shear profile D. The touchdown occurred at 6.0 s at a 
distance of 381.6 m (1,252 ft). The aircraft landed with a pitch attitude 
of 1.727®, a sink rate of 0.658 m/s (2.162 ft/s), and a thrust value of 
10,244 N (2,303 lb). The flare was Initiated at an altitude of 14.78 m 
(48.5 ft). 

Figures 7a to 7c illustrate the performance of the aircraft when It 
flared through the wind profile E. The touchdown occurred at 6.3 s at a 
distance of 399.9 m (1,312 ft). At touchdown the aircraft had a pitch 
attitude of 3.841®, a sink rate of 0.668 m/s (2,194 ft/s), and a thrust 
value of 10,951 N (2,462 lb). The flare was initiated at an altitude of 
15.69 m (51.5 ft). 
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Figures 8 to 10 illustrate the flare performance of the aircraft when 
it encountered a fairly severe shear of +15 kn/30.48 m (100 ft) in the 
horizontal wind velocity. Presently , this is very near the maximum value 
of shear in the horizontal wind component in which the aircraft are allowed 
to attempt landing. 

Figures 8a to 8c show the flare performance of the aircraft when it 
encountered the shear profile F. The touchdowi\ occurred at 5.9 s at a 
distance of 375.5 m (1,232 ft). At touchdown the aircraft had a pitch 
attitude of A. 889*, a sink rate of 0.847 m/s (2.125 ft/o), and a thrust 
value of 10,319 N (2,320 lb). The flare was initiated at a height of 
15.38 m| (50.45 ft), i 

Figures 9a to 9c show the performance in the presence of wind profile 
G. The: touchdown occurred at 5.9 s at a distance of 375.2 m (1,231 ft). 

I 

The aircraft had a pitch attitude of 1.096*, a sink rate of 0.67 m/s (2.201 
ft/s), and a thrust ^value of 11,405 N (2,564 lb) when it touched down. The 
flare had to be initiated at an altitude of 14.49 m (47.55 ft). 

j I 

Figures 10a to, lOc illustrate the flare performance of the aircraft 
when it encountered the shear profile H. In this case, the aircraft 
encountered an increasing headwind of 15 kn/30.48 m (100 ft) and an 
increasing updraft bf 5 kn/30.48 m (100 ft). This was the most severe 
shear profile simulated in the present work. The touchdown occurred at 5.8 
8 at aidistance of 369.4 m (1,212 ft). At touchdown, the aircraft had a 
pitch attitude of 0.2986* (which is barely sufficient to avoid a nose-wheel 
landing), a sink rate of 0.679 m/s (2.227 ft/s), and a thrust value of 
9,857 N (2,216 lb) J Because of the extra lift available from the headwind- 

updraft combination, the aircraft needed to be flared from a much lower 

I i 

attitude of 14.1 m|(46.25 ft). 

l^us, it is seen that even with the use of a single gain setting 
(computed for the shear profile A), the forced-regulator method of 

implementing the control law tends to generate a satisfactory flare 

i i 

trajectory in the presence of widely different wind shear conditions. 

i I 

I 

i CONCLUSIONS 


^rom the timej response curves presented in the previous section for 
the different shear profiles considered, it is clear that the forced- 


I 
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regulator method of ImpleiDontlng the control law developed In an earlier 
section is quite capable of performing the required flare maneuvers in the 
presence of widely differing shear conditions. 

It is noted that, in spite of using a single gain setting (computed 
for the wind shear profile A) for all the cases, three of the touchdown 
conditions, viz the time required for touchdown after initiation of the 
flare, the distance traveled to touchdown, and the thrust value at 
touchdown, were found to be relatively insensitive to the wind conditions. 
Ihe time required to touch down varied from 5. 8 s to 6.3 s. the horizontal 
distance required to Couch down varied from 369.4 m (1,212 ft) to 399.9 m 
(1,312 ft) - a dispersion of +15.2 m (50 ft), which is considered extremely 
good for Che widely differing wind shear conditions studied. The thrust 
value at touchdown varied from a minimum of 9,857 N (2,216 lb) to a maximum 
of 12,557 N (2,823 lb). 

TWo variables, the pitch attitude at touchdown and the altitude at 
which the flare should be initiated [to obtain a sink rate of approximately 
0.64 + 0.04 m/s (2.1 + 0.12 ft/s) at touchdown], were found to be sensitive 
to the wind shear profiles encountered. The value of the pitch attitude at 
touchdown varied from one barely sufficient to a' ''id nosewheel landing 
(0.2896®) to a very robust pitchup (4.889®). It was found that the 
aircraft landed with a lower pitch attitude than the no-wlnd condition 
(fig. 2) for a headwind and a higher attitude for a tailwind. This type of 
behavior is to be expected because of the relative gain (loss) of lift from 
Che resulting headwind (tailwind) for the same pitch attitude. 

The altitude at which Che flare should be initiated Co achieve a sink 
rate of approximately 0.64 m/s (2.1 ft/s) at touchdown also showed some 
variation. The minimum value of this altitude was 14.1 m (46.25 ft) and 
Che maximum value was 15.7 m (51.54 fc). While Chis dispersion in the 
altitude may not be too great, it is considerd Co be a disadvantage. For 
practical purposes, it would be of great convenience to the pilot (or the 
autopilot) Co initiate Che flare maneuver every time from the same decision 
altitude, instead of differing decision altitudes for different wind 
conditions. 

During the course of this project, Che author has reported extensive 
studies of many different types of control laws (refs. I, 2, 8-10, 12). 

The types of control- laws studied varied from time-varying gain to constant 
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gain type, in the absence of any winds to the inclusion of a new wind shear 
model developed by the author (ref. 10), and from the open-loop type of 
control law, and the pure regulator type control law to the forced- 
regulator type of implementation discussed in the present report. ll\e 
forced-regulator method of implementation of the control law was considered 
to have the following advantages: 

(1) ' It was easy to come up with a baseline flare trajectory (in the 
absence of any winds) by imposing two open-loop ramps on the elevator and 
throttle. 

(2) The additional control could be generated by feeding back the 
difference between states due to the winds encountered and the states due 
to no wind conditions. 

(3) Most of the conditions at touchdown (except the pitch and the 
decision height from which flare would be initiated) tere found to be 
insensitive to the widely differing wind shears simulated. This was 
observed in spite of the fact that a single gain-setting was used for all 
the cases. 

The disadvantages of this method of implementation were mainly 
twofold: 

(1) The baseline trajectory in the absence of any winds was found to 
be extremely sensitive to the rate at which the elevator was commanded. A 
very small difference in the command rate resulted in either no touchdown 
(with even a slight rate of Climb achieved) within the duration of time for 
t^hich the simulations were performed or a very hard touchdown, both of 
«rhich are generally unacceptable. 

(2) It is required to select slightly different decision altitudes 
(to initiate the flare maneuver) for different wind profiles. 

The pure regulator type of implementation of the control law was also 
studied extensively by the author (refs. 8-10) and was found to have a 
single, but very important disadvantage - an initial increase in the sink 
rate of the aircraft before it decreased to the desired touchdown value. 

It was determined chat this significant increase in the sink rate was 
solely due Co the way in which the initial conditions were described for 
that problem. The large error in the glide-slope trajectory and the 


18 



I 

I 



shallow trajectory (either triramed or otVierwise) at the desired touchdown 
sinli rate Introduced highly transient behavior in the perturbed trajectory, 
and this was responsible for the significant reversal in the sink rate in 

i 

the Ini^al stage. 

From the extensive study of the different types of control laws fpr 
performing the flare maneuver in the presence of winds, it is evident that 

t 

much work is needed in the area of developing insensitive control laws to 
solve this very difficult problem. Further work is also needed to conduct 
online estimation of the wind parameters and to reduce the sensitivity of 
the touchdown parameters to the differing wind conditions, Ihe results of 
implementing the various types of control laws reported in the present work 
along with that reported in the references should form a base for further 
work in solving this very important problem. 


i 


I 

I 


i 
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Figure 6c. Flare response: wind velocity compo.ients at the wing; wind profile D 
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